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Abstract 

The generahzed Langevin equation is widely used to model the influence of a heat bath upon a 
reactive system. This equation will here be studied from a geometric point of view. A dynamical 
phase space that represents all possible states of the system will be constructed, the generalized 
Langevin equation will be formally rewritten as a pair of coupled ordinary differential equations, 
and the fundamental geometric structures in phase space will be described. It will be shown that 
the phase space itself and its geometric structure depend critically on the preparation of the system: 
A system that is assumed to have been in existence for ever has a larger phase space with a simpler 
structure than a system that is prepared at a finite time. These differences persist even in the 
long-time limit, where one might expect the details of preparation to become irrelevant. 



I. INTRODUCTION 



The influence of a solvent or another complex environment on a chemical motion has 
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routinely been modeled using the Langevin equation or its generalization 
In these equations, complicated many-body effects are approximated by three modifications 
of the equations of motion that describe the intrinsic dynamics of the reactive system: (i) a 
change of the potential energy surface that governs the dynamics, (ii) a stochastic force, and 
(iii) a dissipative friction force. The latter two forces are related by fluctuation-dissipation 
theorems. 

The simplest assumption one can make about the noise is that it is white noise, i.e., the 
strengths of the noise force at different times are statistically independent. This assumption 
requires that the dynamics of the heat bath, which determines the correlation time of the 
fluctuating force, takes place on much shorter time scales than the dynamics of the chosen 
system coordinates, viz., the reaction coordinate and any other solute or solvent coordinate 
coupled to it. In chemical applications, however, the dynamical time scales of the system and 
the bath are usually comparable, so that it is often essential to allow temporal correlations 
of the noise. Fluctuation-dissipation theorems then require that the form of the 

dissipation term must be adapted in such a way that the strength of the damping force 
depends not only on the instantaneous velocity, but on the entire history of the trajectory. 
Thus, the inclusion of correlated noise forces a fundamental change in the mathematical 
structure of the theory: Whereas the Langevin equation with white noise is a stochastic 
differential equation, the generalized Langevin equation that models correlated noise is an 
integro-differential equation that is considerably harder to treat. 

One of the most fundamental properties of a dynamical system is its phase space, which 
represents all possible states of the system in such a way that every phase space point 
corresponds to one and only one state. Recent developments have highlighted the impor- 
tance of adopting a phase space view in dynamical studies of chemical reactivity that are 
based on Transition State Theory (TST): Suitable phase space structures represent the ideal 
recrossing-free dividing surface between reactants and products that was long sought for in 
TST js, 10, 11, 12, 13, 14, 15|. In addition to the dividing surface itself, the phase space view 



of TST provides surfaces in phase space that separate reactive from non-reactive trajectories 

3,y,i2ii. 



and that can be used in an efficient calculation of reaction rates 
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Although all these structures are fundamentally classical objects, they retain their impor 



tance for semiclassical treatments o: 
A recent series of papers 
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29| generalized the phase space view of TST 



to reactions that take place under the influence of time-dependent driving forces or external 
noise. It was shown that all important phase space structures of TST remain intact, but 
move through phase space stochastically. For deterministically driven systems and systems 
under the influence of white noise that are described by the Langevin equation, this work 
could rely on an explicit representation of the phase space: The phase space consists of all 
possible values of coordinates and momenta (or velocities). In these cases, the investigation 
can immediately turn to the problem of describing the dynamical and geometric structures in 
phase space that encode the important features of the dynamics. The generalized Langevin 
equation, however, has a more complicated mathematical structure, and it is not immediately 
obvious what the phase space is or even how many dimensions it has. Several constructions 



of appropriate phase spaces have been suggested in the literature, e.g. in [30 
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will here be shown that the phase space depends critically on the way in which the system is 
prepared. Realistically, every system has come into existence at a certain time in the past, 
but if that time is sufficiently remote, one may simplify the description by assuming that 
the system has existed for ever. Perhaps surprisingly, this assumption changes the phase 
space drastically. A finite time of preparation, no matter how remote, imposes restrictions 
that reduce the dimension of phase space and make its structure more complex. A detailed 
description of how assumptions about preparation infiuence phase space structure even in 
the long term limit, when one might expect them to be forgotten, is the main aim of this 
paper. 

The generalized Langevin equation is usually written in the form 



q = -S/V{q)- [ -fit-s)q{s)ds + ^,{t), 
Jo 



(1) 



where V{q) is a potential energy surface that governs the dynamics of the process under 
study, ^Q,(t) is an external "noise" force exerted by the environment of which typically only 
the statistical properties are known. The integral represents a friction force that the system 
experiences due to its interaction with the environment. The strength of this damping force 
at any given time depends on the prehistory of the system from time t = in a manner 
that is described by the friction kernel 7(t). In thermal equilibrium, the friction kernel is 
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related to the correlation function of the noise by a fluctuation-dissipation theorem 
The form of the friction term used in Eq. assumes that the system is prepared "from 
scratch" at time t = and that the history of the system before that time does not have any 
influence on its future development. That assumption, however, is not always appropriate. 
In some contexts (see, e.g., [25j], where it was convenient to impose a boundary condition at 
t — oo) one might want to assume that the system has been in existence for a long time, 
and thus to replace the GLE ([1]) by 

ij = -VV{q) - f 7(t - s) q{s) ds + ^t). (2) 

J —CXD 

This change might at first glance appear trivial. In particular, one might expect that for 
sufficiently large time t the two versions of the GLE become indistinguishable because the 
difference between the two friction terms will tend to zero. That this conclusion is false 
becomes immediately clear when one considers the number of initial or boundary conditions 
that can be imposed on the solutions of the two GLE: Eq. ([1]) poses an initial- value problem. 
Once the position and velocity at time t = are prescribed, the future of the trajectory 
is completely determined. To determine a trajectory of Eq. ([2]), by contrast, one has to 
prescribe its entire history from t = — cx3 to the present. The prescribed history of the 
trajectory is constrained by the requirement that it must satisfy the GLE. The number of 
free parameters one has to choose can be finite or infinite, but it is always larger than two. 
There is therefore a qualitative difference in dynamics between the two versions of the GLE 
that is manifest in the dimension of their phase spaces and that does not disappear over 
time. 

It is the purpose of the present paper to elucidate the difference in phase space structure 
that is brought about by the seemingly innocuous change in the structure of the damping 
term. The development will be built on a representation of phase space for the GLE that 



was proposed by Martens SSj and that applies to the infinite-time GLE ([2]). A phase space 
for the finite-time GLE ([T]) will be obtained as a time-dependent two-dimensional surface 
in the phase space of the infinite-time GLE, and the dynamics in that surface are given are 
described by time-dependent ordinary differential equations. In this way, it will be shown 
that the dynamics of the infinite-time GLE are simpler than those of its more common 
finite-time counterpart. 

In order to analyze the phase space structures of the GLE, we will for the sake of simplic- 



ity assume that the Brownian particle is moving in a one-dimensional configuration space, 
i.e., the position vector q simplifies to a scalar coordinate q. The formalism to be devel- 
oped generalizes to higher- dimensional problems in a straightforward manner. We will also 
assume that the motion takes place in the vicinity of a parabolic barrier and that therefore 
—W{q) = u^q. This approximation, which covers the case that is most important in reac- 
tion rate theory, allows us to solve the equations of motion explicitly and to check in this 
way that the phase space that is to be constructed has the correct dimension. Higher-order 
corrections to the potential will deform the phase space structures, but leave them quali- 
tatively unchanged, as is discussed in more detail at the end of Section IVIIi Throughout 
most of this paper, we will also leave out the noise terms in Eq. ([T]) and ([2]) and study the 
autonomous equations 

q = ujlq- / -f{t - s) q{s) ds (3) 



^0 

which will henceforth be called the finite-time GLE, and 

q = ujlq- I -f{t - s)q{s)ds, (4) 

J — oo 

which will be referred to as the infinite-time GLE. This omission is justified because we are 
mainly interested in a qualitative theory of the phase space structure. The modifications 
that are required when the noise is included will be described in Section IVIII 

The structure of this paper is as follows: Section [IT] describes the general solution of 
the linearized infinite-time GLE, and thus implicitly its phase space. Secti on [pTI reviews 

1 and that 



the explicit construction of a phase space for the GLE that was given in 
describes the dynamics of the infinite-time GLE. The general solution of the finite-time 
GLE is derived in Section IIVI where it is also shown how the finite-time damping term 
restricts the number of solutions as compared to its infinite-time counterpart. In Section IVT. 
a phase space for the finite-time GLE is constructed as a time-dependent subspace of the 
phase space for the infinite-time GLE. Section IVII illustrates the construction in the simple 
case of an exponentially decaying friction kernel, and it explains why the results obtained 
in this special case are typical of more general friction kernels. Finally, in Section IVIII the 
noise terms is reintroduced into the GLE, and the infiuence of the noise on the phase space 
structure is described. 



5 



II. SOLUTION OF THE INFINITE-TIME GLE 



It is well known (see, e.g., |8|]) that the finite-time GLE ([3]) can be solved by a one-sided 
Laplace transform, the infinite-time GLE (jl]) by a two-sided Laplace transform. However, the 
two-sided Laplace transform is only defined for functions that decay sufficiently fast in both 
the past and the future. In particular, the two-sided Laplace transform of an exponential 
function, which is expected as a solution of the GLE, is undefined. Thus, although the 
calculation formally leads to the correct result, it seems preferable to avoid using Laplace 
transforms. We will here work directly with the GLE. The calculation will here be presented 
in detail in order to highlight the differences between the finite-time and the infinite-time 
GLE. 

If we substitute the ansatz q{t) = e"* into the GLE (jl]), the damping term takes the form 

/7(t — s) q{s) ds = £ 7(t — s) e" ds 
-oo J —oo 

POO 

= ee^* / 7(r) e'^"" dr 







= ee^*7(e), (5) 

where 7(e) is the Laplace transform of the friction kernel 7(t), which is defined by the integral 
in ([5]). The GLE goes over into the nonlinear eigenvalue equation 

^ + e^{e)-ul = Q (6) 

which determines the values of e for which the exponential q{t) = e'^* is a solution of the 
GLE. Because the GLE describes the dynamics in the vicinity of a potential barrier, 
we will in general find one positive real eigenvalue that corresponds to an unstable mode: 
The system slides down the barrier. This eigenvalue is the well-known Grote-Hynes reaction 
frequency 3^ that determines the reaction rate in the space-diffusion limited regime. The 
unstable mode can be interpreted as a collective reaction coordinate that describes correlated 
motion of the reactive subsystem and the bath. It can be constructed either from an explicit 



model of the heat bath 



35 



361] or from a continuum model 37|]. In addition to the unstable 
mode, there will be several eigenvalues with negative real parts that describe stable bath 
modes. Because the eigenvalue equation is analytic and is real on the real axis, complex 
eigenvalues must occur in complex conjugate pairs. 
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The general solution of the infinite-time GLE (jlj) can be written as a linear superposition 

q{t) = Y,o,e^^\ (7) 

i 

where the eigenvalues are the solutions of the eigenvalue equation ([6]) and Cj are arbitrary 
constants. They are subject to the condition that q{t) should be real, i.e., if there are 
complex conjugate eigenvalues, the corresponding prefactors must also be complex conjugate. 
Consequently, the general solution ([7]) has as many free real parameters as there are solutions 
to the eigenvalue equation. This number is the phase space dimension of the infinite-time 
GLE. 

The previous discussion describes the general solution of the GLE ^ correctly if all 
eigenvalues are simple. If the eigenvalue equation ([6]) has multiple solutions, one might 
expect that additional solutions should exist. Indeed, it can easily be checked that the 
ansatz 

q{t) = t^e^* (8) 

with some positive integer k provides a solution of the infinite-time GLE if e is an nth order 
zero of the eigenvalue equation ([6]) and k < n. In this way, an n-fold eigenvalue e gives rise 
to the n solutions 

The total number of independent solutions of the infinite-time GLE (jlj) is therefore the 
number of solutions to the eigenvalue equation (j6j), counted with multiplicities. Thus, this 
number remains unchanged if parameters in the GLE such as the damping strength are 
varied continuously in such a way that eigenvalues coincide. 



III. THE EXTENDED PHASE SPACE 



Martens 



33|, following earlier approaches 



30 



31 



32l |. suggested a construction of the 



dynamical phase space under the assumption that the friction kernel 7(t) satisfies a linear 
differential equation with constant coefficients of the form 
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"+i)(t) + ^a.7«(t)=0. 



(9) 



where 



7«W 



dp 



and ao,...,an are arbitrary real constants. If the friction kernel satisfies this condition, 
one can obtain an explicit representation of the dynamical phase space by introducing the 
customary momentum coordinate p = q and the n + 1 auxiliary variables 

ft 



Co(t) 



t 



l{t - s)p{s) ds, 
7(t — s) p{s) ds, 



(10) 



Cn{t) 



l^''\t - s) p{s) ds. 



These definitions differ from those proposed in [33|] in that the lower limit of integration was 
set to — oo in order to describe the dynamics of the infinite-time GLE. 

With the help of the variables d the GLE can be rewritten as a system of linear ordinary 
differential equations 



q = p, 



p = ujlq- Co, 

(. = 7«(0)p + 0+i forO<z<n, 

n 



(11) 



i=0 



In the last of these equations, the differential equation ([9]) of the friction kernel has been 
used. If these coordinates are gathered into the vector 

(q\ 



P 
Co 



(12) 



the equation system ( jTTj) takes the form 



Mz 



(13) 
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with the matrix 



M 



/o 


1 








.. 









-1 





.. 


. 





7(0) 





1 


.. 


. 





7(0) 








1 . . 


. 



(14) 



\q yn)(0) -ao -a, -as ... -a^/ 
The differential equations ( fTTi) or ( fT3l) do not exphcitly contain memory effects. The future of 
a trajectory is uniquely determined once the present values of the coordinates g,p, Co) • ■ • , Cn 
have been prescribed. The auxiliary coordinates C,i encode all the relevant information 
about the history of a trajectory. In order to obtain a closed system of equations for finitely 
many coordinates, it is essential that the friction kernel satisfies the condition ([9]). If no 
such condition holds, one obtains instead a differential equation system that involves all 
coordinates Q for < i < oo. The phase space of the system is then infinite dimensional. 
Even in this case, however, the construction shows that the phase space dimension is always 
countably infinite. 

The formulation (fTTj) of the GLE in terms of a system of linear differential equations with 
constant coefficients can be solved with the usual techniques. A solution 



Z{t) = ZqC 



et 



with exponential time dependence exists if e is an eigenvalue of the coefficient matrix M. 
Martens 33|] has shown that these eigenvalues are exactly the solutions of Eq. ([6]), so that 
his construction of the phase space, where it is applicable, leads to the same solutions as the 
explicit calculation of Sec. [Tll 



IV. SOLUTION OF THE FINITE-TIME GLE 



The finite-time GLE (|3]) poses and initial-value problem: A trajectory is completely 
determined once the position g(0) and velocity q{0) have been prescribed at time t = 0. 
This observation indicates already that the general solution of the finite-time GLE cannot 
be of the form (^^, which contains too many free parameters. Nevertheless, one would expect 
that the eigenvalues obtained from Eq. ([6]) should still determine the dynamics. Indeed, it 
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will be shown in the following that the phase space of the finite-time GLE can be regarded as 
a time-dependent surface of the extended phase space that describes the infinite-time GLE. 
It will also be demonstrated how the extended phase space must be restricted to model the 
dynamics of the finite-time GLE. 

If we once again make an exponential ansatz q{t) = e*^* and substitute into the finite-time 
GLE, we find that the friction term takes the form 



In comparison to the earlier result for the infinite-time GLE, there is a correction term 
that arises from the finite lower limit of integration. Even if e is chosen to be a solution to 
the nonlinear eigenvalue equation ([6]), the exponential ansatz will fail to provide a solution 
to the GLE, unless the boundary term vanishes for all times t, which is impossible. The 
finite-time GLE therefore does not have solutions with simple exponential time dependence. 
If more general solutions of the form ([71) 



This equation, which must be satisfied for all t, provides constraints on the parameters Cj. 
It therefore reduces the dimension of the phase space below the number (including multi- 
plicities) of eigenvalues ej. 

The discussion of the previous section allows us to rewrite the boundary condition f|T5|) 
in a more suggestive manner. Solutions of the finite-time GLE are, by assumption, defined 
only for times t > 0. However, a solution of the known analytical form ([7]) can easily be 
extended to negative times. If we take the liberty to do this, (ITSl) takes the form 




j 



are permitted, the condition that the boundary term should vanish reads 




(15) 




(16) 
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which should again hold for all times t. For a trajectory that satisfies this condition, the 
finite-time and infinite-time friction terms yield the same result. This condition therefore 
demands that the history of the trajectory before time t = should have no infiuence on the 
friction force at any time t > 0, which is precisely the assumption by which the finite-time 
GLE differs from the infinte-time GLE. 

By setting t = in Eq. (fT6|) . we obtain the condition Co(0) = 0, with the auxiliary 
variable Co introduced in the previous section. If 7(t) is sufficiently differentiable, we can 
take derivatives of the condition (fT6!) to obtain 

0(0) = for all i > 0. (17) 

If the friction kernel satisfies a differential equation of the form (Q, only the first r?, + 1 of 
the conditions (II 7p are independent, and the validity of the remaining conditions is enforced 
by the differential equation. 

This observation allows us to interpret the physical origin of the boundary condition (ITSjl : 
The auxiliary coordinates Q encode all the information about the past of a trajectory that 
is necessary to predict its future. A point (g,p, Co, • • • i Cn) in the extended phase space 
corresponds to a possible state of the system if and only if a past trajectory can be found 
for which all variables take the prescribed values. At first glance, it might appear obvious 
that this should always be possible, at least for the infinite-time GLE, because one has 
the freedom to adjust the (infinitely many) values of the function q{s) for s from — oo to 
the present, and one has to match only a finite or countably infinite number of prescribed 
values for the phase space coordinates. However, if the function q{s) is to describe a possible 
history of the system, it must satisfy the GLE for all past instances of time. This condition 
drastically reduces the freedom in the choice of a past trajectory, and it is not obvious what 
precisely the consequences of this requirement will be. 

In the case of the infinite-time GLE, the differential equations (ITTl) allow one to calculate 
the past as well as the future of a trajectory for an arbitrary initial condition {q,p, Co) ■ • • i Cn) 
in the extended phase space. This observation shows that for every such initial condition it 
is always possible to identify a suitable past trajectory q{t), and this trajectory is unique. 
The extended phase space therefore is the dynamical phase space of the system in which 
every point corresponds to one and only one possible trajectory. 

The situation is drastically different in the case of the finite-time GLE. While it is clearly 
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possible to prescribe the values g(0) and p{0) at the initial time t = 0, one cannot choose 
arbitrary values for the d at this time because the trajectory does not have a past that could 
be adjusted to the prescribed Q. Indeed, from the modified definition (which is the original 



definition of Ref. 



m 



m= / i^'\t - s) p{s) ds (18) 

that is appropriate to the finite-time GLE, it is obvious that 0(0) must always be zero and 
cannot be chosen arbitrarily. We are thus led back to the condition (|T7]) . (Note that the 
difference between the original definition (|TOl) and the modified definition (|T8l) is immaterial 
under the constraint (fT6l)). As was to be expected, the only initial conditions that can be 
imposed for the finite-time GLE are the values q{0) and p{0). 



V. THE PHASE SPACE OF THE FINITE-TIME GLE 

The extended phase space described in Section [TTll cannot be regarded as the true dy- 
namical phase space of the finite-time GLE because its dimension is too high. Instead, the 
true phase space is a two-dimensional surface within the extended phase space that is char- 
acterized by the constraints f|T71) . An important subtlety to note about these constraints is 
that they impose conditions on the initial values of the auxiliary coordinates Q, not on the 
current values. To decide if a point in the extended phase space is in the true phase space 
at time t, one has to follow the trajectory through this point back to time t = and then 
check the initial values Ci{0). For a given point the initial values of will depend on t. As 
a consequence, the true phase space is a time-dependent surface S(t) of the extended phase 
space. This surface, which for the linear GLE ([3]) is a plane, will in the following be referred 
to as the phase plane, as opposed to the extended phase space in which it is embedded. 
The observation that the finite-time GLE does not permit solutions with purely exponential 
time dependence is intimately related to this time dependence of the phase plane. 

We will now derive equations of motion for the dynamics in the phase plane. It will 
be shown that these equations are still linear, but explicitly time dependent. Consider the 
dynamics in the extended phase space that is described by Eq. f fT3|) 

z = Mz 

with the matrix M of Eq. (fT4|) . (For the sake of definiteness, we will focus on the case 
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in which the extended phase space has finite dimension. Similar arguments apply to the 
infinite-dimensional case.) Let U{t) be the time evolution operator of Eq. f[T3|) . i.e., let the 
dynamics be given by 

z{t) = U{t)z{0). (19) 

Let the operator P(t) project a point {q,p, Co, • • • ) in the extended phase space onto the 
initial values (Co(0), Ci(0), • • • ) of the auxiliary coordinates. Points in the phase plane E(t) 
that represents the dynamical phase space at time t are then described by the condition 

zej: ^ P{t)z = 0. (20) 

At time t = 0, this projector has the form 

^0 1 ...^ 
F(0) = 1... 

V • • • 

Because the P{t) projects onto initial values, it satisfies 

P(0)z(0) = P{t)z{t) = P{t)U{t)z{0) 

or 

P(0) = P{t)U{t) = const. (21) 
Taking the time derivative of this equation, we find 

P{t) = -P{t)U{t)U-\t) = -P{t)M 

because U = MU. 

To introduce a coordinate system in the phase plane S, we choose two orthogonal unit 
vectors Vi{t) and V2{t) in S. These vectors have to be time-dependent in such a way that 
they stay in S at all times, but this requirement does not determine the time-dependence 
completely. We can impose the condition that their derivative Vi be perpendicular to the 
plane S, i.e. 

Vi-Vj = forz,j = l,2. (22) 
This choice also ensures that the orthonormality conditions 

vl = vl = l, vi-V2 = (23) 
13 



are satisfied at all times if they are satisfied at t = 0. 

Because the vectors Vi are to be in E, they must satisfy the constraint 



P{t)vi{t) = 0. 



Taking the time derivative of this equation, we obtain 



P{t)vi{t) = Pit) Mviit), 



whence it follows that 



Vi = MVi + \iiVi + Xi2V2 



with scalar functions Xijif). The orthogonality conditions (l22l) determine these functions as 



These equations of motion determine the coordinate vectors Vi{t) completely once i'i(O) 
and 1^2(0) have been arbitrarily chosen as an orthonormal system in S(0). Note that the 
right-hand side of Eq. (1241) is simply the component of Mvi that is perpendicular to 
Thus, the coordinate vectors are subject to the dynamics f[T^ . only corrected to satisfy the 
constraint (1221) . 

Once the coordinate system in the phase plane has been chosen, an arbitrary point in 
the plane can be described with two coordinates y\ and 2/2 by 



and the dynamics in the phase plane will be described by equations of motion for the yi. 
These equations can be obtained from f[T^ by noting that 

yi = Vi - Mz 



Xij = Vj ■ Mvi 



which gives 



Vi = MVi - {Vi ■ MVi)vi - {V2 ■ MVi)v2. 



(24) 



z = yiVi + y2V2 



(25) 



{vi ■ Mvi)yi + {vi ■ Mv2)y2 



(26) 



or 




(27) 
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where 

(M),., J'-'^'-'C"-*^-'! (28) 

y{V2-MVi) {V2-MV2) 

is the projection of the matrix M that describes the dynamics in the extended phase space 
onto the instantaneous position of the phase plane E. This projection is exphcitly time de- 
pendent, and it must be because the plane S is itself time-dependent. As a consequence, the 
equations of motion (I27p do not permit solutions with purely exponential time dependence, 
although they are linear. This confirms the observation made in Section IIVI The devel- 
opment of the present section gives a geometric interpretation of this observation, and it 
provides an explicit construction of the phase plane that represents the true two-dimensional 
dynamical phase space of the finite-time GLE. 



VI. EXAMPLE: EXPONENTIAL MEMORY FRICTION 

As an example, we will study the case of an exponentially decaying memory kernel 

7(t) = ^e-*/- (29) 

T 

with a characteristic memory time r and a damping strength 7. This friction kernel is 
normalized such that 



7(t) dt = 7. 

-00 

It satisfies a first-order differential equation of the form iQ and therefore leads to a three- 
dimensional extended phase space, the smallest dimension possible. 



331] . There is always 



The geometry of the extended phase space was described in detail in 
one positive real eigenvalue. The two remaining eigenvalues have negative real parts and 
correspond to damped bath modes. They are either both real or form a complex conjugate 
pair. Excitations in the bath modes either decay monotonically or perform damped oscilla- 
tions, respectively. The parameter regions that lead to each type of behavior are illustrated 
in Fig. [H 

The dynamics in the true phase space is obtained by projecting from the extended phase 
space onto the phase plane The motion of this plane is itself determined by the 

dynamics in the extended phase space. For this reason, in the long time limit the phase plane 
will, loosely speaking, be aligned with the unstable and the slowest-decaying stable mode. In 
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FIG. 1: Regions in the parameter space of the infinite-time GLE with exponential damping (|29p 
that lead to monotonic decay or damped oscillations in the bath modes. 

1.0 \^ 
0.5 ^ 




FIG. 2: Eigenvalues of the projected dynamics described by {M)Y:(t) ^ function of time for 
7 = O.StJb s^nd uj]jT = 0.2. Dotted lines indicate the two largest eigenvalues of the dynamics in the 
extended phase space. 

the case of three real eigenvalues in the extended phase space, this is literally true: For long 
times, one of the coordinate vectors Vi is parallel to the unstable eigenvector, the second 
is parallel to the eigenvector with largest negative eigenvalue, subject to the condition that 
the two vectors must be perpendicular. The influence of the smallest eigenvalue disappears 
in this limit. That this is actually the case is shown by the numerical example in Fig. El 
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The figure shows the eigenvalues of the matrix {M)Y,(t) that describes the instantaneous 
dynamics in the phase plane as a function of time. At t = 0, these eigenvalues are itUb- 
This must be the case for any choice of the friction kernel 7(t), as can easily be verified from 
Eq. f[T^ or from the fact that the finite-time damping term in Eq. Q is zero for t = 0. 

The observation that both eigenvalues and eigenvectors of the dynamics within S(t) 
converge to fixed limits makes it easy to describe the long-time dynamics. There will be a 
stable manifold, i.e. a set of trajectories that approach the origin as t ^ oo. At sufficiently 
large t, this stable manifold will coincide with the subspace spanned by the stable eigenvector 
within the phase plane. However, the picture is more complicated at early times, when the 
time-dependent eigenvectors have not yet reached their limiting values. The stable manifold 
is itself time- dependent, and at early times it will not be aligned with the instantaneous 
stable eigenvector. Instead, at these times a full solution of the time-dependent equations of 
motion is necessary to identify those trajectories that will ultimately approach the origin. 

Although there is an unstable eigenvector within the phase plane, an unstable manifold 
can, strictly speaking, not be defined. Such a manifold is generally defined as the set of 
trajectories that approach the origin as t — — oo, but this limit cannot be taken for the 
finite-time GLE. Nevertheless, it is clear that for t ^ oo all trajectories outside the stable 
manifold will exhibit exponential instability along the direction of the unstable eigenvector, 
so that the space spanned by this eigenvector plays the role of an unstable manifold for 
practical purposes. 

A more complex situation, illustrated in Fig. [3l arises if the stable eigenvalues of the 
extended dynamics form a complex conjugate pair. As before, at time t = the eigen- 
values of the dynamics within the phase plane take the values ±0;^. For large time, the 
unstable eigenvalue approaches the unstable eigenvalue of the extended dynamics and one 
of the coordinate vectors Vi aligns with the corresponding unstable eigenvector. The second 
eigenvalue, however, does not approach a fixed limit, but a limit cycle. Because the two sta- 
ble eigenvalues in the extended phase space describe oscillatory dynamics, the phase plane 
is forced to rotate around the unstable direction. The projected matrix (M)s(i) takes 
the same value twice in every rotation period because the orientation of is irrelevant. 
Indeed, it can be checked numerically that the oscillation period in Fig. [3] is half the period 
described by the imaginary part of the complex eigenvalues of the extended dynamics. 

In the case of three real eigenvalues, it is easy to conclude that for large times the 
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FIG. 3: Eigenvalues of the projected dynamics described by {M)Y,{t) as a function of time for 
7 = O.Swb and = 2. Dotted line indicates the unstable eigenvalues of the dynamics in the 
extended phase space. The other two eigenvalues form a complex conjugate pair. 

phase plane must be spanned by the two eigenvectors that correspond to the leading 
eigenvalues. No such simple prediction is possible in the case of complex eigenvalues. It is 
clear that the phase plane will contain the unstable eigenvector and will carry out a rotation 
in the plane spanned by the two stable eigenvectors. The frequency of that rotation is 
given by the imaginary part of the eigenvalues, but its phase can only be obtained from 
an explicit solution of the equations of motion, not from qualitative considerations. The 
coordinate direction that is aligned with the unstable eigenvector will still play the role of 
a (quasi-) unstable manifold. The instantaneous stable eigenvector of the dynamics within 
the phase plane, however, will be periodically time-dependent even in the long-time limit, 
and it is therefore not obvious where the stable manifold of the origin will be. Again, this 
can, of course, be determined from a full solution of the equations of motion, but not from 
a simple qualitative argument. 

The exponentially decaying friction kernel that was discussed in the previous section 
represents the simplest possible case because it leads to the lowest possible dimension of the 
extended phase space. However, the characteristics of the dynamics within the phase plane 
that were found in this special case are typical of the general case because any fast decaying 
additional bath modes that are present in the extended phase space lose their influence on 
the dynamics in the long-time limit, as was shown in the first scenario of Sec. IVII 

The motion of the phase plane, which is represented by the motion of the coordinate 
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vectors Vi, is determined by the dynamics in the extended phase space. In the long-time 
hmit, one of the coordinate vectors will always be aligned with the unstable eigenvector. The 
second coordinate vector will be determined by the eigenvector or eigenvectors corresponding 
to the bath modes that decay most slowly. If there is a single such eigenvector with a negative 
real eigenvalue, the coordinate vector will be aligned with it (as far as possible subject to the 
constraint that it must be perpendicular to the first coordinate vector), and the dynamics 
in the phase plane will qualitatively show the behavior of the first scenario in Sec. I VII If 
the slowest-decaying mode is a damped oscillation, the phase plane will rotate around the 
unstable eigenvector as described by this mode, and the dynamics within the phase plane 
will follow the pattern of the second scenario in Sec. IVI[ 

The two cases discussed above will therefore correctly describe the dynamics in the phase 
plane of the finite-time GLE in the long-time limit for arbitrary friction kernel, except pos- 
sibly in degenerate cases where several decaying or oscillatory bath modes have eigenvalues 
with the same real part. For short times, of course, the influence of the strongly damped 
modes is not negligible. If one's goal is, for example, to find stable manifolds at time t = 0, 
it is necessary to take full account of all phase space dimensions, and the simplifications that 
are afforded by the long-time limit are not available. 

VII. INCLUSION OF THE NOISE 

The most important feature of the generalized Langevin equation ([1]) and ([2]) is the noise 
term. It might appear odd, therefore, that throughout this work the noise term has been 
neglected. It is crucial to see how the inclusion of the noise will change the picture that has 
so far been presented. 

It is straightforward to incorporate the noise into the phase space formulation of the 
infinite-time GLE. A noise term will appear in the second equation of the system ( fTTil . 
which will then read 

P = ^Iq- Co + ^a{t). 

The equations of motion (fT3l) are replaced by 

z = Mz + Ut), (30) 
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where 



( \ 





V / 

The noise changes individual trajectories, but the dynamics still take place in the extended 
phase space spanned by Co, ■ ■ ■ , Cn, or possibly q,-p and infinitely many Ci. 

The case of the finite-time GLE is, once again, not as simple. A trajectory is determined 
by two initial conditions and the realization of the noise that is acting upon it. The phase 
space of the finite-time GLE must therefore still be two-dimensional, but if it is regarded as 
embedded into the extended phase space, its embedding will depend on the noise. 

Specifically, the autonomous time evolution (IT^ in the extended phase space is replaced 
by a time evolution 

z(t) = U^(t-zm (31) 

that depends explicitly on the realization a of the noise. The time evolution is no longer 
linear, and therefore cannot be represented by a matrix, because the noisy equations of 
motion are inhomogeneous. As in the noiseless situation, the phase space of the finite-time 
GLE at time t = is the plane S(0) that contains all points in the extended phase space 
that satisfy Ci = 0. At later times, the phase space contains all points that evolve out these 
initial conditions, i.e. the phase space of the finite-time GLE is obtained through the time 
evolution of S(0) as 

= /7,(t;S(0)). (32) 

Significantly, the phase plane that is obtained in this way is not only time-dependent, as it 
was in the noiseless version of the finite-time GLE, it also depends on the realization a of 
the noise. 

The nonlinear time evolution fl3Tl) for the noisy dynamics on a harmonic barrier can be 



obtained explicitly by the same approach that was used in [24, |25| to obtain a geometric 
description of Transition State Theory for the Langevin equation: The solution of the noisy 
equation of motion (15U]1 is split into two components 



(33) 
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Here zl^{t) is a particular solution of the noisy equation (l30l) that remains in the vicinity of 
the barrier top at all times and is called the Transition State trajectory. As the subscript a 
indicates, it depends on the noise. Explicit formulas for zl^{t) as an integral over the noise 
are given in 2J, |25| (where the formulas require knowledge of the noise at times t < 0, it 
can be chosen arbitrarily for the finite-time case). The relative coordinate Az{t) satisfies 
the noiseless equation of motion 

Az = M Az 
and therefore follows the noiseless time evolution 

Az{t) = U{t) Az(0) 

with the same linear time evolution operator as in Eq. fll9l) . We can then write 

z{t) = Uit)AziO) + zUt) 

= uit) (z(o)-4(o)) + 4(t) 

= U^{t;z{0)), 

which gives an explicit expression for the nonlinear time-evolution operator defined by 
Eq. ( 13T1) . It can be solved for 

z{0) = U-\t)z{t)-U-\t)zi{t) + zUO). 

The condition for a point z{t) to be in the phase plane of the finite-time GLE at time t can 
then be stated as 

= P(0)z(0) 
= P(t)2(t)-P(t)4(t)+P(0)4(0), 
= P{t)z{t)-s^{t), 

where the time evolution Eq. fl2Tl) of the projector P{t) has been used. This condition, which 
characterizes the phase plane Sq,(^), can be restated as 

Pit) z{t) = s^{t) (34) 

with the vector 

s4t) = P(t)zi(t)-P(0)4(0). (35) 
21 



It describes a time-dependent plane that is parallel to the phase plane T,(t) for the noiseless 
case that was defined in fl20l) . but it is shifted away from the origin by an amount that is 
stochastically time dependent. As anticipated, therefore, the phase space of the finite-time 
GLE for any given realization of the noise is a two-dimensional plane within the extended 
phase space. 

An arbitrary vector within the phase plane can be written as 

z = yivi + y2V2 + Sa (36) 

where the two vectors Vi(t) and V2{t) span the noiseless phase plane S(t) (see Section |V|) 
and the reference vector Sa{t) satisfies P{t)Sa{t) = Sa{t). Eq. fl30l) then translates into the 
equations of motion 

yi = Vi- {Mv^)y, + ■ {Mv2)y2 + ■ ^^{t) - v., ■ S^{t) (37) 

or 

= {Mhit) r']+ {Ut))m - {Sa{t))m, (38) 
) \y2 J 

where (M)s(f) is, as in (!28|) . the projection of the matrix M onto the plane S(t) and the 
last two terms denote the projections of the given vectors onto that plane. In comparison to 
the noiseless equation of motion ( 1271) . the term {$,a{'t))'E{t) describes precisely the impact of 
the noise on the dynamics that one would expect. The unexpected term —{Sa{t))Y:{t) arises 
from a motion of the reference vector Sa{t) along the instantaneous direction of the phase 
plane. This term can be eliminated by a suitable redefinition of Sa{t)'- If an arbitrary vector 
in S(t) is added to Sa{t), the shifted phase plane SQ,(t) remains unchanged. 

The formalism presented in this section can directly be generalized to systems with an- 
harmonic barriers. In this case, the finite-time GLE ([T]) will lift to a nonlinear equation of 
motion 

z = f{z) (39) 

in the extended phase space. As before, this equation of motion must be solved subject 
to the constraint Ci(O) = on the auxiliary variables, and this constraint defines a two- 
dimensional plain S(0) that represents the phase space at time t = 0. The solution of the 
equation of motion (l39i) is described by a nonlinear and stochastic time evolution operator 
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z{t) = Ua{t; z{0)) as in Eq. (13T]) . The phase space SQ,(t) at a time t > is given by the 
time evolution of the initial phase plane S(0) according to Eq. fl32l) 

Because the time evolution operator is nonlinear, the phase space T,a{t) will in general not 
be a plane. For an arbitrary anharmonic potential in the GLE, it will usually be impossible 
to obtain explicit formulas for either the time evolution operator Ua or the phase space 
T,a{t). Nevertheless, in qualitative terms the situation remains unchanged: The dynamical 
phase space of the finite-time GLE is represented by a two-dimensional submanifold of the 
extended phase space that is moving stochastically through the larger extended phase space. 

VIII. CONCLUDING REMARKS 

The development presented here demonstrates that there are important differences be- 
tween the dynamics of the finite-time and infinite-time GLE. While the infinite-time GLE 
allows for trajectories with purely exponential time dependence, as one would expect for a 
linear equation of motion, the finite-time GLE imposes boundary conditions that prohibit 
such behavior. These conditions reduce the phase space dimension for the finite-time GLE 
to two. From a geometric point of view, the dynamics of the finite-time and infinite-time 
GLE differ drastically. The infinite-time GLE can be rewritten as a set of time-independent 
ordinary differential equations in a phase space with dimension between three and infinity. 
The finite-time GLE, by contrast, has a two-dimensional phase space, the dynamics in which 
is described by differential equations that are time-dependent in a complicated way. 

In spite of these dramatic differences in the geometric framework, the dynamics of the 
finite-time and infinite-time GLE become similar when one studies individual trajectories in 
the long-time limit. In both cases, the behavior of a trajectory in this limit is determined 
by one unstable eigenvalue (which is the same in the extended and two-dimensional phase 
spaces). As the boundary condition imposed by the finite-time GLE recedes into the past, its 
influence on a particular trajectory becomes less and less palpable, although the differences 
between the respective phase spaces, which represent the totality of possible trajectories. 
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remain for arbitrarily long times. 
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